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NATO  UNCLASSIFIED 


AN  UNRESTRICTED  LINEAR  RANDOM  WALK  WITH  NEGATIVE 
EXPONENTIALLY  DISTRIBUTED  STEP  LENGTHS 


By 


B.W.  Conolly 


ABSTRACT 


An  account  is  given  of  the  theory  of  a  doubly  infinite  linear  random  walk  in 
which  step  lengths  have  a  negative  exponential  distribution  and  the  direction  of 
each  step  is  not  necessarily  equiprobable.  The  problem  of  first  passage  time 
is  also  studied.  The  theory  was  developed  in  connection  with  a  study  of  random 
linear  anti-submarine  patrols. 


NATO  UNCLASSIFIED 
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INTRODUCTION 


The  unrestricted  linear  random  walk  described  in  text  books  is  usually  a  study 
of  the  motion  of  an  imaginary  particle  that  roams  backwards  and  forwards  along 
the  doubly  infinite  x-axis  according  to  definite  probability  laws.  The  axis  is 
divided  into  equal  intervals  of  length  h,  say,  and  the  particle  proceeds  in  leaps 
and  bounds,  the  lengths  of  which  are  integral  multiples  of  h.  At  the  end  of  a 
leap,  which,  henceforward,  to  conform  with  convention,  we  shall  call  a  'step'  , 
the  particle  chooses  the  length  and  direction  of  its  next  step.  The  probability 
distribution  associated  with  this  choice  is  prescribed  in  advance  and  remains 
the  same  throughout  the  walk. 

The  most  usual  topics  of  investigation  are,  for  a  given  starting  point: 

th 

a.  the  probability  that  at  the  end  of  the  n  step  the  particle  is  at  the 
point  with  x -coordinate  mh,  where  m  is  a  positive  or  negative  integer; 

th 

b.  the  probability  that  at  the  end  of,  or  during,  the  n  step  the  particle 
reaches  its  starting  point  for  the  first  time  since  the  walk  began.  (This  is  the 
so-called  'first  passage  problem'.') 

If  the  particle  moves  with  constant  speed,  steps  will  occupy  varying  times 
proportional  to  their  lengths.  On  the  other  hand,  if  it  is  supposed  that  a  step 
occupies  constant  time,  the  particle  has  to  move  with  variable  velocity.  In 
either  case  a  plot  of  the  particle's  position  as  a  function  of  time  produces  a  sort 
of  random  noise  record  consisting  of  a  sequence  of  straight  line  segments,  and 
in  this  connection  topic  (b)  above  has  an  obvious  relation  to  the  problem  of  the 
interval  of  time  between  successive  zeros  of  a  random  function. 
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The  random  walk  described  and  investigated  in  this  report  arose  from  the  study  of 
a  family  of  tactical  problems.  It  differs  from  the  discrete  step  random  walk  in 
that  the  lengths  of  steps  are  continuous.  At  the  end  of  a  step  the  particle  is  given 
a  choice  as  to  whether  to  move  to  the  right  with  probability  p,  or  to  the  left  with 
probability  (=  1-p),  and  the  length  of  the  next  step  is  governed  by  the  negative 
exponential  distribution.  Thus,  the  probability  that  a  step  has  length  between 
s  and  s+ds  is  ywe  ds. 

The  tactical  application  of  the  first  passage  time  aspect  of  this  walk  is  described, 
together  with  the  skeleton  of  the  derivation  of  the  relevant  probabilities,  in  the 
classified  Ref.  1.  It  is  felt  desirable  to  set  down  an  unclassified  account  of  the 
derivation  of  the  formulae,  not  only  for  the  first  passage  time  problem,  but 
also  for  the  problem  of  where  the  particle  is  at  the  end  of  a  given  step.  A 
similar,  but  in  some  aspects  abbreviated,  account  has  also  been  published  by 
the  author  in  Ref.  2. 
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1. 


POSITION  OF  THE  PARTICLE  AT  THE  END  OF  THE 


th 


n 


STEP 


A  particle  is  considered  that  executes  a  random  walk  on  a  line  in  such  a  way 

S 

that  the  probability  that  a  step  has  a  length  between  s  and  s+ds  is  yue  ^  ds. 
When  the  particle  reaches  the  end  of  a  step,  it  has  probability  p  of  moving  to 
the  right  at  the  next  step,  and  q=l-p  of  moving  to  the  left. 


We  study  here  the  joint  probability  and  probability  densities  y) 

L  (x,  y)  that  the  particle,  which  begins  its  walk  at  the  point  X  with  coordinate 
^  th 

X,  finds  itself  at  the  end  of  the  n  step  at  a  point  Y  in  the  small  interval 
(y,  y+dy).  When  y  x,  y)  used  to  denote  this  probability;  when  y— x, 

L^(x,  y)  is  the  corresponding  form. 


The  fundamental  equations  satisfied  by  y)  and  L^(x,  y)  are  as  follows; 


Rn^x,  y)  =  qyu 


Jo 


y+s)  e  '^^ds  +  p/U 


/^y-x 

7 

Jo 


s)  e  ds  + 


y-s)  e  ds 


(Eq.  1) 


L^(x,  y)  =  pytA  J*  L^_^(x,  y-s)  e  '^^ds  ^  j  L^_^(x,  yfs)  e  '^^ds  + 


+ 


I  /  R  j(x,  y+s)  e  ^^ds 

Jx-J 


(Eq.  2) 


The  first  term  in  Eq.  1  deals  with  the  case  when  the  particle  is  to  the  right 
of  Y  at  the  end  of  the  (n-1)^^  step.  The  second  and  third  terms  arise  when  the 
particle  is  to  the  left  of  Y  at  the  end  of  the  (n-1)^^  step,  the  integral  from  0 
to  y-x  dealing  with  the  case  when  it  is  to  the  right  of  X,  and  the  final  term  when 
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th 

the  end  of  the  (n-1)  step  is  to  the  left  of  X.  A  similar  explanation  can  be 


given  for  Eq.  2. 

We  can  formulate  similar  equations  by  considering  the  particle's  first  step. 
Suppose  this  is  to  the  left  (probability  q),  and  to  the  point  within  the  small 
interval  (x-s,  x-s  — ds).  Then  the  particle  must,  during  the  next  n-1  steps, 
march  from  x-s  to  y.  The  resulting  integral  is: 


If  the  first  step  is  to  the  right,  two  cases  arise  according  to  whether  it  is  to 
the  right  of  Y  or  to  the  left.  Thus,  altogether 


oo 


r 


R  (x, 
n 


y)  =  q/^f 

•/© 


e  R^_j(x+s,  y)  ds  + 


(Eq.  3) 


A  comparison  between  Eqs.  1  &  3  shows  that  R^(x.  y)  is  a  function  of  y-x. 
Similarly,  L,^(x,  y)  is  a  function  of  x-y. 


An  evaluation  by  first  principles  of  Rj(x,  y),  Lj(x,  y),etc,  gives: 


(Eq.  4) 


/ 


(Eq.  5) 
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^2^^’  ^  P  e  [  q  +  P/A  (y  -  x)]  ,  (Eq.  6) 

^2^^’  q/A-e^^^^^  j^p+q/A  (x-y)  J  ,  (Eq.  7) 


and  these  lead  us  to  consider  the  functions 


^n^^"  y);  defined  by 


R  (x,  y)  =  e  ^  (y-x) 

n  '  n 


(Eq.  8) 


L  (x, 
n 


y) 


■/<  (x-y) 


^  (x-y) 
n 


(Eq.  9) 


Evidently,  either  ^  ^ 

other  by  interchanging  p  and 


n 

q, 


(and  R 

n 

X  and 


and  L  )  can  be  transformed  into  the 
n 

y- 


Finally  we  introduce  the  generating  functions 


/O  {  S  ,  t)  = 


A(  5  .  t)  = 


and  the  Laplace  transforms 


r  (z)  = 
n 


(  5  )  dS 

n 


(Eq.  10) 


(Eq.  11) 


(Eq.  12) 


6 


L 


.  oo 

I  (-z)  -  !  e 

"  ■  0 


X  (  6  )  d6 

n 


r(z,  t)  = 


CO 


-z5 


(  S  ,  t)  dS  , 


iiz,  t)  =  X  (  5  .  t)  d^ 

Jq 


Clearly 


r(z,  t)  = 


r  (z)  t” 
n  / 


nil 


with  a  similar  definition  for  tiz,  t). 


(Eq.  13) 


(Eq.  14) 


(Eq.  15) 


(Eq.  16) 


Since  at  the  end  of  each  step  the  particle  must  be  somewhere,  we  have  for 
each  n  ; 


oo 

R^(x,  y)  dy  + 


L^(x.  y)  dy  =  1. 


(Eq;  17) 


or 


or 


X  iS  )  dS 

n 


=  1, 


(Eq.  18) 


r  {  M  )  +  i  {  M  )  =  1 
n  '  n  ' 


(Eq.  19) 
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Hence,  for  |  t  |  <  1  and  Rl  z  ~  ytA.  ,  r(z,  t)  +  €  (z,  t)  is  analytic,  and,  in 
particular , 


r(  /<  ,  i)  +  t  { ,  t)  =  t/(l  -  t)  • 


(Eq,  20) 


The  integral  equation  for  ^  6  )  deriving  from  Eq.  1  is 

/>n<  ^  )  =  1/  /  /J^_j(6  +  s)ds+p^  / 

Jq  J.Q 


-s)ds  + 


+  pyu  e 


2/a6 


/ 


oo 


(Eq.  21) 


■2^ 


X  ,  (s  -  S  )  ds  . 
n-1 


Applying  the  Laplace  transformation,  we  obtain 


r  (z)  -  qyU 


n 


and 


fr  (z)-r  .(2/^  r  1 

Ljlzl - fr  ,(z)  +  l  ,(2/.)] 


(Eq.  22) 


(2^  -  z) 


(Eq.  23) 


These  are  simple  difference  equations,  and  it  is  easy  to  show,  using  Eq.  19,  that 


r  (z)  =  /A. 
n  ' 


-z) 


(Eq.  24) 
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r^(2^  -z) 


(Eq.  25) 


^  (z) 
n 


n 


The  second  term  on  the  right-hand  sides  of  Eqs.  24  &  25  are  particular  solutions 
of  Eqs.  22  &  23,  which  have  to  be  included  to  complete  the  general  solution. 
Since  1*^(2 yw.  )  is  certainly  not  infinite  (see  above)  one  must  expect  that 
t^2 -z)  contains  terms  that  annihilate  the  apparent  poles,  and  it  seems  a 
reasonable  conjecture  that  both  and  ^^(z)  are  n^^  degree  polynomials 

in  descending  powers  of  z. 


n 


Multiplying  Eqs.  22  &  25  by  t  ,  adding,  and  using  the  results 


r  j(z) 


we  obtain 

qyK  tzr  (2yw. .  t)  -pyUi  t(2yp>.  -z)/ (2^  ,  t) -p^  t(2yu  -z) 

r(z,t)  =  - - - .  (Eq.  26) 

2  2 
z  z(2-qt  +  pt)  +  2p yA  t 


with  a  similar  formula  for  •£.(z,  t),  in  which  p  and  q  are  interchanged.  We 
denote  the  zeros  of  the  denominator  of  Eq.  26  by  and  Then 


"^1 

=  /*! 

[  1  +  1  (p-q) 

“2 

'A  1 

[  1  +  1  (p-q)  t-i^ 

“  1 

+  1  (q-p)  t+i^ 

“2 

1^1  +  1  (q-p)  t-]^ 

y 


(Eq.  27) 
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where 


=  (1  -  |t)^  -  pqt^  . 


(Eq.  28) 


Now  and  are  not  less  than  yu  for  all  p,  and  t<  1.  Hence,  since 

r(z,  t)  and  /(z,  t)  are  analytic  for  R1  z  and  t  <  1,  we  must  arrange 

for  the  numerator  of  Eq.  26  and  its  companion  expression  for  i,  (z,  t)  to  vanish 

n 


respectively  at  z  =  and  z  = 


Thus,  noting  that 


Y  °  “'2 

Y  -"'i  '“2  • 


(Eq.  29) 


we  obtain 

q  coj  r(2^  ,  t)  "  pw'2  ^(2/K.  ,  t) 
qco^  r(2yix  ,  t)  -  pw'^  •£(2^  ,  t) 


-q 

2  j 


and  hence 

r(2yu  ,  t) 


pa'l  +  q 

q(co'i  coj  -o'j  coj) 


(Eq.  30) 


(Eq.  31) 


^(2^  ,  t) 


P co'  +  q  CO.  co„ 
^2  2^1  2 


P(“’l  “1  '“’2  “2^ 


(Eq.  32) 
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Substituting  in  Eq.  26  we  obtain 


r(z,  t) 


t(p  w'j  + 
2R 


1 

(z  -U2) 


^(z,  t) 


t(qco^  +  p  w'^) 
2R 


These  may  be  inverted  easily  and  give 


/O  (  &  .t)  = 


t<p“'i  +  q"2>  6«2 


2R 


A  (  6  ,t)  =  - - 


t(qco  +  p  w’  ) 


2'  6w'< 


2R 


and^  returning  to  the  original  generating  functions^ 


R(x,  y,  t) 


Y. 

nil 


t(pco'^  +  qw 
2R 


-(y-x)(/A 

e 


L(x,  y,  t) 


t(qw^  +  pw'2) 
2R 


-(x-y)(/A 

e 


It  is  of  particular  interest  to  study  the  densities  y),  E^(x, 

n  and  various  y.  The  simplest  case  is  that  in  which  p  =  q  =  ^ 


(Eq.  33) 

(Eq.  34) 


(Eq.  35) 

(Eq.  36) 


(Eq.  37) 


(Eq.  38) 

y)  for  fixed 
,  when  the 
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particle  is  as  likely  to  move  to  the  right  as  to  the  left  during  each  step. 
R(x,  y,  t)  and  L(x,  y  t)  are  then  identical,  and,  from  Eqs.  37,  27  &  28, 


R(x,  y,  t) 


(Eq.  39) 


where  5  -  y-x.  Now  R  (x,  y),  being  the  coefficient  of  t^  in  R(x,  y,  t),  is 

n 

given  by 


R  (x,  y) 

ll 


_1 

n! 


d"  R 


dt 


n 


_ 1_ 

2  TT  V 


I  R(x,  y,  z) 

0 


(Eq.  40) 


where  C  is  a  simple  closed  contour  surrounding  the  origin,  such  that  R(x,  y,  z) 
is  analytic  in  z,  inside  and  on  C. 


We  make  the  transformation 


Then 


R  (x,  y) 
n 


2  TT  4. 


-/A.S  (1 

e 


-5  ) 


dl 

J^(2-5)”  ’ 


(Eq.  41) 
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where  C’  is  a  simple  closed  contour  containing  ^  =  0  and  excluding  S  =  2. 
Evaluating  the  residue  at  J*  =0  we  have 


V-)  =  I 

2  ^ 


ji"CL 

e  V  i2^S  f 


rl 


/  2n-r-2\ 

V  ^-1  A 


(Eq.  42) 


where  6  =  y“X. 


The  remaining  paragraphs  are  intended  to  provide  a  check  on  Eq.  42. 


Since 


f. 


oo 


r 


y)  dy  +  /  L^(x,  y)  dy  =  1 

X  Aoo 

for  every  n,  and  since  in  this  case  (p  =  q)  L^(x,  y)  and  y)  are  identical, 

except  that  in  the  case  of  the  former,  6  in  Eq.  42  stands  for  x-y,  we  would 
expect  each  of  the  two  terms  on  the  left-hand  side  to  be  equal  to  one  half.  From 
Eq.  42, 


OO 


n— 1 


R  (x,  y)  dy  =  - 
'X  2 


1 _  ^  /2n-r-2\  _ 1 

n  Vn-1  y  2^-1 

TmO 


-r 


(Eq.  43) 


If  we  write 


oo 


(i-i  )(2-  :J)^ 


-I 

riiO 


k  <■’’  t 

r 


(Eq.  44) 


then  it  is  easy  to  confirm  that 


r 


R_^(x,  y)  dy  = 


(n) 


(Eq.  45) 
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We  note  the  identity 


(1-?  ){2-S) 


n  +  1 


(Eq.  46) 


Then 


(n) 

^n-1 


2n-l  ^  1  ^  ,  (n+1) 

n  I  „ 2n  n-1 


(Eq.  47) 


We  now  proceed  to  show  that  the  right-hand  side  of  Eq.  47  is  equal  to  k 


(n+1) 


n 


From  Eq.  44 

,  (n+1) 

n-1 

Then 


.n+1 


i  (”■•) 

BmO 


1  /2n-lV|^(n.. 

^  2n  ^  n  y  n-1 


(n+1) 


1 _  /  2n  \  (n+1) 

2n-l(nj  n-1  ^ 

1 _  ^  /2n-s\  1 

,n.l  4  I  n  y  ^n-s 


2 

(n+1) 


n 


i.  e.  k 


(n) 

n-1 


(n+1) 


n 


(Eq.  48) 
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This  common  value  is,  in  particular,  true  for  n  =  1,  and  evaluation  of 
as  the  coefficient  °  in  [  (1-  r )  (  2  -  r )  shows  the  common 

value  to  be  1/2.  Thus 


\(x,  y) 


dy  = 


as  expected,  and  a  check  is  provided  on  Eq.  42. 
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2. 


DISTANCE  TRAVELLED  DURING  FIRST  PASSAGE  TO  ORIGIN 


In  this  section  we  consider  the  problem  of  the  distance  travelled  by  the  particle 
from  its  initial  point  X  (coordinate  x  ^  0),  until  it  reaches  the  origin  for  the 
first  time.  To  be  more  precise,  we  consider  the  probability  y)  dy 

(y  >  X  ^  0)  that  simultaneously: 

(a)  the  particle  reaches  or  crosses  the  origin  for  the  first  time  at  the 

th 

end  of,  or  during  its  n  step; 

(b)  the  total  distance  traversed  by  the  particle  up  to  the  moment  when 
it  first  reaches  the  origin  lies  between  y  and  y  +  dy. 

It  is  convenient  to  introduce  two  new  probabilities:  y)  dy  and  y)  ^y* 

The  probability  y)  dy  has  a  definition  identical  with  that  of  y)  dy» 

except  that  the  particle's  first  step  is  to  the  right.  The  probability  y)  dy 

corresponds  to  a  first  step  to  the  left.  Clearly  yh  y^> 

are  all  zero  for  y  ^  x  and  n  =  1,  while, for  y^  x  and  n  ^  2  , 

y)  =  p  f^(x,  y)  +  q  g^(x,  y)  •  (Eq.  49) 

The  case  y  =  x  is  a  special  one.  The  particle  can  return  to  the  origin  in  n 

steps  covering  a  total  distance  x  only  if  it  moves  always  to  the  left.  We 
denote  the  corresponding  probability  by  v^(x),  and  then 


Vj(x)  =  q  e  ^  ^  , 


(Eq.  50) 
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and  generally 


£ 


V  (x)  =  q  /  V  ,  (x  -  G)  dG 

n  /  n-1 


(Eq.  51) 


which  gives 


V  (x)  =  q  e 

n 


■/*  X  (q/A  x) 


n-1 


(n  -  1  )  I 


(Eq.  52) 


The  generating  function  (g.f.  )  of  v^(x)  is  given  by 


.n 


v(x,  t)  =  'y  t“ 


I 

n^l 


(Eq.  53) 
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3. 


FORMULATION  OF  EQUATIONS 


We  consider  first  the  case  x  >  0  and  write  down  integro -difference  equations 
for  y)  and  g^(x,  y),  for  y  ^  x.  We  have 


1 

tn<x,y)=j/*e 


'n-1 

J  0 


-vuG 

e  '  ^(x+G,^  y  -  G)  dG  • 


(Eq.  54) 


The  first  term  in  Eq.  54  arises  because  the  particle's  farthest  possible  first 

step  to  the  right  cannot  carry  it  beyond  the  small  interval  at  the  point  with 

coordinate  “  (x  +  y)  (a  step  length  of  r  (y  -  x)),  and  a  direct  passage  to  or 

beyond  the  origin  in  n-1  steps  must  follow.  The  integral  term  arises  to  cover 

all  the  possible  steps  to  a  point  in  the  interval  (x  +  G,  x  +  G  +  dG),  where 

0  ^  G  ^  ^  iy  -  x).  Thence  the  particle  starts  another  walk,  the  first  step  of 

which  can  be  to  the  right  or  the  left,  and  it  must  reach  or  overstep  the  origin 

th 

for  the  first  time  at  or  during  the  (n-1)  step,  covering  a  total  distance  to  the 
origin  of  an  amount  between  y-G  and  y  -  G  +  dy.  Similarly 


e  ^  u^  ^(x  -  G,  y-G)  dG 


(Eq.  55) 


By  substituting  Eqs.  54  &  55  in  Eq.  49  we  obtain  a  homogeneous  integro - 
difference  equation  for  y^- 
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An  obvious  way  of  converting  integro -difference  equations  into  integral  equations 
is  to  introduce  g.fs.  Thus,  let 


F(x.  y.  t) 


G(x,  y,  t) 


U(x,  y,  t) 


E 


ni2 


E 


ni:2 


E 


n^2 


.n 


t  f^(x,  y) 


n 


t  y) 


n 


t  u^(x,  y) 


/ 


The  analogues  to  Eqs.  49,  54,  &  55  are  then 

U(x,  y,  t)  =  p  F(x,  y,  t)  +  q  G(x,  y,  t)  » 
1  .  ^  /*-J-(y-x) 


1  —/^(y-x) 

F(x,  y,  t)  =  -/At  e  Vj^^x+y),  t 


J ] 

J  0 


(Eq.  56) 


(Eq.  57) 


e"'^*^U(x+G,  y-G,  t)  dG, 


£ 


(Eq.  58) 

X 

-/aG 


G(x,  y,  t)  =  /.At  j  e  ^  U(x-G,  y-G  ,  t)  dG 

'O 


(Eq.  59) 


In  passing,  it  is  worth  mentioning  that  F(x,  y,  t)  and  G(x,  y,  t)  satisfy  the 
simultaneous  differential  equations 


Qf 

3x 

3f 

ay 

II 

1 

-/A  t  U, 

(Eq. 

60) 

dx 

^G 

dy 

+  /A.G  = 

+/A.t  U, 

(Eq. 

61) 
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which  imply  that  all  three  of  F,  G,  and  U  satisfy  the  telegrapher’s  equation 


F  =  0  , 

where 

a'  =  ^  +/«t(p  -  q)  -/<(2  -t)  -  y  (1  - 1) 

Ox  Oy 

(Eq.  62) 

and  the  omitted  arguments  x,  y,  and  t  are  understood. 


We  now  proceed  to  determine  F,  G  and  U  from  Eqs.  57,  58,  &  59.  We  first 
introduce  the  function  F  (x,  z,  t),  G(x,  z,  t),  U  (x,  z,  t),  and  V  (z,  t),  where 


F  (x,  z,  t)  = 


/ 


oo 


-zy 


F(x,  y,  t)  dy 


(Eq.  63) 


G  and  U  are  similarly  defined,  and 


V(z,  t)  = 


L 


QO 

V(y,  t)  dy 


qt  ^-x(a-q^t) 


a 


(Eq.  64) 


We  then  observe  that 


f. 


e  dy  I  e  U(x-G,  y-G,  t)dG  =  e  /  e^^  U(s,  z,  t)  ds,  (Eq.  65) 


.oo 

e  dy  /  e'^^  U(x+G,  y-G,t)dG  =  e®^/  e^®  U  (s,  z,  t)  ds. 


r 


(Eq.  66) 
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where 


a  =  /M-  +  z  • 


(Eq.  67) 


Then 


ax  — 


Fix,  z,  t)  =  e  V  (a+z,  t)  +  yu.  t 


U(s,  z,  t)  ds,  (Eq.  68) 


/ 

0 


_  #  3.  S 

G(x.  z.  t)  =  ^  te  /  e  U  (s,  z,  t)  ds. 


(Eq.  69) 


whence,  using  the  transformed  form  of  Eq.  57 


ax  — 


U  (x,  z,  t)  =  ^  pte  V  (a+z,  t)  +  ^  pte 


£ 


oo 

ax  f  -as  -  xNj  , 

e  U (s,  z,  t)ds  + 


+  ju  qte  e^^  U  (s,  z,  t)  ds  , 

Jo 


(Eq.  70) 


To  solve  Eq.  70  for  U  we  introduce  the  further  Laplace  transform 


oo 

U  (co,  z,  t)  =  /  e  ^^U(x,  z,  t)  dx. 


I 


(Eq.  71) 


which,  applied  to  Eq.  70,  leads  to 


|■qt(a)-a)+U(a,  z,  t)(a-^^  qt)(a+a)-^  qt) 


r  2  2  1.  1 

Tco  +yutco(p-q) -(a  -ayut)^  (a yu.qt)  (a +cj  qt) 


(Eq.  72) 
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Now  u  (x,  y)— >-0  as  x-*»  oo  for  any  n,  since  it  is  clearly  impossible  for  the 
n 

particle  to  reach  the  origin  in  a  finite  number  of  steps  from  an  infinite  distance. 
It  follows  that  U  (w,  z,  t)  is  an  analytic  function  of  w,  at  least  for  Rico  ^  0  and 

0  <  t  <  1. 

The  poles  of  the  right  hand  side  of  Eq.  72  are  at 

w  =  yuqt  -  a,  (Eq.  73) 


0) 


0)  -  0)^ 
2 


-  j  ^t(p  -  q)  +  t^{p  -  q)^  +  4 

J  F/*  t  (p  -  q)  -  (p  -  q)^  +  4 


(Eq.  74) 


(Eq.  75) 


where 

^2  2 

Tv  =  a  -  ayw  t  .  (Eq.  76) 

Taking  t  ^  1  we  find  that 

R1  (  yu  qt  -  a)  <  0, 

so  that  the  pole  co  =Qa.  qt  -  ^is  admissible.  Also  R1  co^  ^  0  for  all  p  and  q 
while  Rj  co^  y  0.  Hence  the  numerator  of  the  square  bracket  in  Eq.  72  must 
vanish  for  co  =  co^.  This  fact  determines  U  (a,  z,  t)  and  we  obtain 
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U  (co,  z,  t)  = 


yupqt  (a  +  cj) 


(co  -  (a  +  cl)^  -/*Aqt)  (a  +  co  -^qt) 


/*-pqt 


(a+o)  -/Aqt)(a  +  co  -^qt)  c 

X  ^ 


cOj  +  a 


CO  -  CO,  CO  - 


q/^ 


t  +  a 


(Eq.  77) 


Thus,  on  inversion,  we  have 


U (x,  z,  t)  = 


r  (CO 

L  1 


COj^X 


+  a)  e  -  ^qte 


-(a -/A.  qt)x 


]■ 


(Eq.  78) 


and,  upon  substitution  in  Eqs.  68  &  69, 


a  +  CO  \  CO  X 

F  (x,  z,  t)  =  I  -  1  t  e 


a  -  CO 


1 


(Eq.  79) 


c^-.  X  ,  .. 

.  z,  t)  =  t  e  -  t  e'^^ 


G(x 


(Eq.  80) 


Since  the  Laplace  transform  U*(x,  z,  t)  of  U(x,  y,  t)  is  given  by 


-zx 


U*(x,  z,  t)  =  e  V(x,  t)  +  U  (x,  z,  t)  , 


(Eq.  81) 


we  have 


(a  +  CO  )  CO,  X 

(x,  z,  t)  =  -  e 


(Eq.  82) 
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Similarly,  the  Laplace  transforms  of  F(x,  y,  t)  and  G(x,  y,  t)  are  given  by 


F*(x,  z,  t) 


(11^ 
V  ^  -  “l 


-) 


t  e 


COj^X 


(Eq.  83) 


GHx,  z,  t) 


t  e 


COj^X 


(Eq.  84) 


Recovery  of  F(x,  y,  t),  G(x,  y,  t),  and  U(x,  y,  t)  can  now  be  carried  out  by 
inverting  the  transforms.  Thus  we  obtain  for  y  ^  x  ^  0 


G(x,  y,  t)  = 


P-q/*. 


xt 


-/*y(i-|t)-|/it(p-q)x 

e 


(Eq.  85) 


and 


U(x,  y,  t)=  e 


(y  +  x) 


V 


pq(y^-x^) 


+ 


(Eq.  86) 


F(x,  y,  t)  can  obviously  be  obtained  immediately  from  U  and  G.  If  desired, 
the  densities  u^(x,  y),  etc.  can  be  extracted  from  the  expansion  of  Eq.  86  in 
powers  of  t. 
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In  the  case  where  the  particle  starts  from  the  origin  (x  =  0),  the  problem  of  first 
passage  back  to  the  origin  concerns  a  path  in  which  the  first  n-1  steps  are  either 
entirely  to  the  right  of  the  origin,  or  entirely  to  the  left.  In  the  former  case,  and 
with  obvious  notation. 


1 

•j/*y 


y-G)  dG  . 


(Eq.  87) 


Corresponding  to  a  walk  to  the  left  of  the  origin  there  is  a  probability  density 

f  '  (0,  y)  similar  in  all  respects  to  f  (0,  y),  except  that  p  and  q  are 
n  n 

interchanged.  We  have 


u  (0,  y)  =  pf  (0,  y)  +  qf  (0,  y)  • 
n  n  *^  n 


(Eq.  88) 


The  same  kind  of  treatment  as  for  x  ^  0  yields  the  Laplace  transform 
U  (0,  z,  t)  of  the  generating  function  U(0,  x,  t).  We  obtain 


U  (0,  z,  t)  =  (  2  a  -/At  -  s)  , 


(Eq.  89) 


where 


2 

s 


2  2,  ,2  .2 

/At  (p-q)  +4^ 


2 

)\  =  a  -  a^  t  . 

a  =  /A  +  z  # 
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On  inversion  this  gives 


U(0,  y.  t) 


-  2  -  t) 

e  Ij  (/^ty  1^)  . 


(Eq.  90) 
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